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Extensive numerical calculations show that the capacitance of back-gated nanowires with various 
degrees of dielectric embeddings is accurately described with an effective dielectric constant as long 
as the difference between the dielectric thickness and the gate-nanowire distance is held constant. 
This is valid for dielectrics with permittivities ranging from simple air to water. However, due to 
screening the scaling is not valid if the dielectric lies down well below the nanowire. Moreover, when 
only the dielectric thickness varies the capacitance characteristics are S-shaped with three distinct 
regions, of which only the first two can be used for dielectric sensing. The first region is almost 
linear while the middle region, with a span of two diameters around the center of the nanowire, is 
the most sensitive. 
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I. INTRODUCTION 

The continuous progress encountered in the interdis¬ 
ciplinary field of nanotechnology has made possible the 
fabrication, characterization, and utilization of a plethora 
of nanosized structures in 1, 2, and 3 dimensions. One 
of the most studied categories of nanostructures is that 
of nanowires (NWs) and, in particular, semiconductor 
NWs, in which the charged carriers are confined in 2 
dimensions and can eventually move freely in the third 
direction. The semiconductor NWs offer complex build¬ 
ing blocks for devices and structures with functionalities 
used in fields like nanoelectronics or photovoltaics [T]. 
Their specific shape makes them ideally as active elec¬ 
tronic interfaces with biological cells; hence the active 
structures with a NW field effect transistor (FET) are of 
great promise not only as electronic devices but also as 
interfaces and sensors in biological systems [2]. In ad¬ 
dition, electrical measurements and probing of biological 
cells and tissues with the help of NW-FET display high 
signal-to-noise ratios mostly because of protrusion capa¬ 
bilities and small active areas of those devices [3] . 

A metallic-like behavior in NWs is obtained with a 
doping as large as carriers/cm^ where the De¬ 

bye screening length is of a few nanometers (i. e. much 
smaller than the NW diameter) [4]. In this case the 
metallic behavior of the NW brings the NW-FET in 
the linear regime where there is a linear relationship be¬ 
tween the conductance and the gate voltage. In the linear 
regime there is a definite relation between the back-gate 
capacitance Cq of the NW-FET and the effective mobil¬ 
ity which is one of the key parameters characterizing 
the charge transport in the FETs. Thus if we assume that 
the NW-FET is of length Lq^ under a drain-source volt¬ 
age VdSi a gate-source voltage Vqs the source-drain 
current Ids' is given by the following expression [5] 


Ids = RfeCg{Vgs - yt)VDslL%, ( 1 ) 

where Vt is the threshold voltage that is defined as 


the gate voltage that fully depletes the free carriers in 
the NW. The effective mobility iipE can be calculated 
from the slope of the current versus gate voltage or the 
transconductance Qm = dlos/QV gs with the following 
formula 


Ree 



CgVds 


( 2 ) 


Equation ([^ is used for estimating the charge carrier 
mobility in NWs and shows the key role of the back-gate 
capacitance. Most of the estimations of the back-gate 
capacitance are based on the assumption that the NW is 
a long cylinder at the distance t above a metallic plane 
associated with the back gate. The capacitance of such 
a back-gated system is given by 


_ 27rgQgy.I/G 

^ cosh“^(t/R) 

where R is the radius of the cylinder, and Sr is the relative 
dielectric constant of the dielectric in which the system 
is supposed to be totally embedded. 



FIG. 1: Cross-section view of the NW. The NW is depicted 
as a circle of radius R, G is the gate, t is the distance from the 
NW center to the gate and h is the thickness of the dielectric 
represented as dotted lines. 

Typically, the configuration of a NW-FET is shown in 
Fig. I^a), where the NW is separated from the gate G by 
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a dielectric slab of thickness h = t — R. The capacitance 
formula given by Eq. has been commonly used in the 
evaluation of back-gate capacitance of NW-FETs [ 6 ]-[ 8 ]. 
On the other hand, Eq. © is correct if only the entire 
space is a dielectric of permittivity 5 ^. As a result it was 
noticed in several works that, for real devices, formula 
given by Eq. is valid if the relative permittivity of the 
dielectric slab Sr is replaced by an effective permittivity 
Seff that is approximately half (more precisely 58%) of 
Sr when the dielectric is silicone dioxide [iiini. 

In the present work we perform an extensive study 
by considering not only different dielectric slab thick¬ 
nesses or embeddings but also different dielectric con¬ 
stants. Some of the present results have been previously 
reported in a recent conference paper m- Two of those 
embeddings are depicted in Fig. [^b) (the NW is more 
than half-buried by the dielectri^ and Fig. [^c) with 
the NW fully buried just below the dielectric-air inter¬ 
face. We will show that Eq. © may be used to describe 
the back-gate capacitance of the NW for arbitrary slab 
embeddings. This description of the capacitance is deter¬ 
mined hy t — h that is the difference between the distance 
NW-gate and the dielectric thickness. Thus, if t — h is 
constant, Eq. ® is valid with an appropriate Seff that 
depends not only on t — h but also on Sr- Deviation from 
this behavior occurs when the dielectric slab is thinner, 
i. e., (t — h)/i? ^ 1. We further analyze capacitance 
changing with respect to dielectric thickness in order to 
assess the sensing capabilities of liquid levels in microflu¬ 
idic systems. We have found that with respect to the 
dielectric thickness the capacitance characteristics have 
an S-like shape with three distinct regions. The capaci¬ 
tance varies almost linearly in the first region, slowly in 
the third one, and rapidly in the middle region that is 
located around the NW. 

The paper has the following structure. In the second 
section we present the problem and the method of reso¬ 
lution. The third section is dedicated to the main results 
and the last section will summarize the conclusions. 


II. PRELIMINARIES AND THE METHOD 

Our problem resides in solving the Laplace equation 
with Dirichlet boundary conditions, i. e., a fixed poten¬ 
tial V on the surface of cylinder and with a grounded 
back-gate. The capacitance can be calculated numer¬ 
ically by a variety of methods including finite eiement 
methods [T 2 j or boundary eiement methods [13]. The 
boundary eiement method is the finite eiement version of 
the boundary integrai method which comes from poten- 
tiai theory [T4|. In the static limit metallic regions are 
considered to totally screen the electric fields thus those 
regions are equipotential regions. This is not valid when, 
instead of static fields, the electromagnetic fields are con¬ 
sidered, such that metals have a finite dieiectric permit¬ 
tivity at opticai frequencies. Thus, for metaiiic nanos¬ 
tructures with features much smaiier than the waveiength 


of the incoming eiectromagnetic wave the response of a 
metaiiic system is given by the same Lapiace equation 
but with different boundary conditions (basicaiiy there 
are Neumann boundary conditions) [T5l[T6|. The method 
is a spectrai approach to the boundary integrai equation 
that turns out to obtain simuitaneousiy aiso the eiectro- 
static capacitance of metaiiic nanostructures HZ]. 

Capacitance and the eiectrostatic response of finite 
cyiinders have been studied in HU and in m , respec- 
tiveiy, where explicit expressions for both capacitance 
and electric polarizability of finite cyiinders have been 
found. On the other hand, an infinite iong cyiinder with 
a grounded pianar back-gate and homogeneousiy embed¬ 
ded in a dieiectric of reiative permittivity Sr has aiso an 
explicit capacitance that is given by Eq. One cal¬ 

culation procedure is based on the bipolar coordinates 
miiin]. In the usual setup with the NW perpendicular 
on the (x,y) plane, the bipolar coordinates given by 


_ asinh(r]) 

cosh(?7)—cos(^) 
_ Qsin(^) 

^ cosh(?7)—cos(^) 

Z = ^. 


(4) 


generate the equipotential surfaces and the field lines of 
a cylinder charged under a potential V with respect to 
the grounded planar gate HB. The equation of the NW 
surface is 77 = 770 , while the equation of the planar gate is 
77 = 0. We notice that the ratio R/t from Fig. [^is t/R = 
cosh(77o) In the (x,y) plane the new coordinates ( 77 , 
are shown in Fig. Thus, the surfaces 77 = constant are 
cylinders that surround our NW m- They turn out to be 
equipotential surfaces of a charged cylinder including also 
the surface 77 = 0 that is our planar gate. The surfaces ^ = 
constant are also cylinders and define the fieid iines that 
start up normaiiy on the NW and end up aiso normaiiy 
on the pianar gate m- 



FIG. 2 : The bipolar coordinates ( 77 , ^). Surfaces 77 = constant 
surround the NW and ^ = constant start up on NW and end 
up on the back-gate G. 

We have shown that in the case of electrostatic capac¬ 
itance once the equipotential surfaces and field lines are 
known one can estimate straightforward the capacitance 
of the system HZ]. Accordingly, the capacitance of the 
back-gated NW is 
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where and hz are the Lame coefficients of the 

transformation @ 

h = _a_ 

^V cosh(?7)—cos(^) 

cosh(?7)—cos(^) 

hz = 1 . 

Combining ^ with ^ will lead to 


^ _ 27rSoSrLG 
Vo 


( 7 ) 


which is just Eq. ^ since 770 = cosh“^(t/i?). Similar for¬ 
mulae can be obtained for two non-concentric cylinders 
or for two parallel cylinders HD. 




FIG. 3: 2D cross-section of equipotential surfaces calculated 
(a) with a finite box and (b) with infinite boundary elements. 
The dielectric is shown by the horizontal line intersecting the 
NW which is placed in the middle of the pictures. 

Slight modifications of the shapes of nanostructures 
induce minute changes on capacitance m or on higher 


excitations modes [20l[2T]. On the other hand, the back- 
gated NWs are not infinitely long and have a certain de¬ 
gree of doping, effects that cannot be analytically quan¬ 
tified. Extensive numerical calculations have shown that 
with a doping as much as 5x10^^ cm“^ and with an as¬ 
pect ratio Lg/R as large as 100 the behavior of the semi¬ 
conductor NW is that of infinitely long metallic NW [22] . 
Eurthermore, the capacitance of the back-gated NW with 
a finite but otherwise arbitrary thickness of the dielec¬ 
tric cannot have a readily analytic formula like that of 
the homogeneous case. Hence we invoke a fully numer¬ 
ical procedure to calculate the NW capacitance in vari¬ 
ous dielectric embeddings. Our calculations are based on 
ANSYS which is a finite element based multiphysics soft¬ 
ware program [23|. In ANSYS there are several methods 
for computing the electrostatic capacitance of a metallic 
system. The first and the simplest method is a pure fi¬ 
nite element method in a finite computation box, hence 
the equipotential surfaces are enforced to close on the 
boundaries of the computation box (Eig. [3}i). Another 
method is the method that uses infinite boundary ele¬ 
ments (Eig. [§). As we can see from the figure the 
infinite boundary elements ensure a more physical ap¬ 
pearance of the equipotential surfaces with respect to a 
finite computation box by considering the same computer 
overhead. Therefore, the second method is more accurate 
and faster in convergence than the first one. Neverthe¬ 
less we have used the first method but we have tuned it 
with respect to the second one by making the results of 
the two methods be apart by less than 0.5%. The major 
reason for this choice was the fact that the first method 
is more suitable to be used in ANSYS scripts of repeated 
calculations. 


III. RESULTS 


Our main results are presented in Eigs. i E and 
m The dielectric constants considered in this work are 
Sr = 3.9 (i. e. Si02), = 10(a lower limit of high-/^ di¬ 

electrics), and Sr = SI (the dielectric constant of water). 
In Eig. I^we present the NW capacitance as a function of 
t/R with t/R >5 (the typical setting for a NW-EET) and 
with different dielectric embeddings. The level of dielec¬ 
tric embedding was established in such a way that t — h 
is constant. In this setting the capacitance of the NWs 
shows good scaling properties according to Eq. ^ with 
an effective dielectric permittivity Seff for each value of 
t — h. The particular case {t — h)/R=l, in which the di¬ 
electric separates the NW and the gate plane, is plotted 
with upward oriented triangles and is also fitted with an 
effective dielectric permittivity. The results of fitting are 
plotted with solid lines showing a quite good superposi¬ 
tion on the calculated curves. The case {t — h)/R=l is the 
standard configuration of the NW-EET and was studied 
in |9] where an effective permittivity Seff = 0.57esiO2 
was given. We have obtained a value Seff = 0.bSesiO2' 
In addition, when the high-/i: dielectric Hf02was used 
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FIG. 4: NW capacitance versus tjR when t — h is kept 
constant. The curve /i=0 is the case of air as a dielectric, 
while the curve h =inhnity is the case of completely dielectric 
Filings. The cases (a), (b), and (c) depicted in Fig. are, 
respectively, (t — h)/R =1, -0.2, -1. The ht of curves for 
(t — =1 is plotted with solid lines. 


{Sr = 25) the effective permittivity was about 34% of 
that of the dielectric [9]. Our calculations indicate that 
the effective permittivities are 43 percent and 5.5 per¬ 
cent of permittivities of the dielectric with Sr = 10 and 
Sr = 81, respectively. One can further notice that for 
{t — h)/R varying between -1 and 1 the Ceff varies non¬ 
linear ly with respect to the variation of Sr- 


The scaling is still valid for dielectric permittivities as 
large as that of water in spite of the fact that at quite 
large dielectric constant the dielectrics should behave as 
a metal. Nevertheless the metallic behavior (manifested 
as strong screening) of the dielectric sets in when (t — 
h)/R 1 {i. e., the NW is well above the dielectric 
layer) and the scaling is no longer valid. One can also 
notice that at {t — h)/R = 2.4 for both Cr = 10 and 
Sr = SI the capacitance is almost constant, a feature of 
metallic behavior. Fig. illustrates more clearly that 
only for (t — h)/R > 1 (Fig. If) the dielectric screens 
the back gate hence the scaling described above cannot 
longer be valid. In addition, for dielectric embeddings 
that completely cover the NW the capacitance goes very 
slowly to the value of complete filling of the space with 
dielectric (Fig. [^). 

We will further analyze and discuss the previous scal¬ 
ing behavior for certain applications used to determine 
the dielectric thickness by electrical measurements using 
NW-FET. A potential application one can think of is the 
estimation of the liquid height in microfluidic systems. 
To keep the discussion as simple as possible we suppose 
that the surface of the NW is passivated and the liquid is 
totally wetting the NW surface. In Fig. [^we plotted the 
NW capacitance with respect to thickness h at various 
values of t/R. In fact the curves plotted in Fig. [^give 
information about Seff discussed previously. The capac¬ 
itance is ^'-shaped with three distinct regions. For small 
thicknesses the capacitance varies linearly as it was ex¬ 
pected on perturbative grounds. The size of this region 
increases but the slope of the curves decrease with the 
increase of t/R. Hence, at t/R = 40 the size of this lin¬ 
ear region is significantly larger as one can see it from 
Fig. § In this region, in order to exploit its linearity 
for sensing, one has to trade off between the size of the 
linear region and the slope of capacitance with respect 
to dielectric thickness h. By inspecting Fig. it seems 
that t/R = 10 — 11 ensures an optimum sensing in this 
region. The second region displays the largest capaci¬ 
tance variation with greater slopes for greater dielectric 
permittivities. The region spans approximately four NW 
radiuses around the NW center and can be used as a 
proximity sensor in which the capacitance increases con¬ 
siderably when the liquid levels are in the vicinity of the 
NW. Finally, the third region of deeper dielectric embed¬ 
ding the capacitance varies smoothly and saturates to its 
asymptotic value of total embedding. 


IV. CONCLUSIONS 

In the present work we have studied the back-gate ca¬ 
pacitance variation of a nanowire with respect to various 
levels of dielectric embedding. The standard configura¬ 
tion in the field effect transistor setting is that in which 
the dielectric separates the nanowire from the gate plane. 
It was previously shown that this standard geometric ar¬ 
rangement has a scaling behavior straightforwardly con- 
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FIG. 5: Cross-section of the equipotential surfaces for dielec¬ 
tric: (a) below the NW, (b) partially burying the NW, and 
(c) above the NW. Only (a) is able to screen the back gate. 
We have used Sr = 10. 


nected to the back-gated nanowire embedded in a homo¬ 
geneous dielectric. The scaling is made with an effective 
dielectric constant that strongly depends on the permit¬ 
tivity of the dielectric. In our paper we have shown nu¬ 
merically that this scaling is valid for various degrees of 
dielectric embedding granted the fact that the difference 
between the dielectric thickness and the gate-nanowire 
distance is constant. The scaling is dependent on the 
difference between the dielectric thickness and the gate- 
nanowire distance and on the permittivity of the dielec¬ 
tric. However, the scaling is not valid for dielectric thick- 



0l I . I t , , , . 

2 4 6 8 10 12 14 16 18 

h/R 





2 4 6 8 10 12 14 16 18 

h/R 


FIG. 6: NW capacitance versus h/R when t/R is: 5, 6, 7, 8, 
11, and 40. The corresponding vertical lines indicate the ratio 
t/R in order to delimitate the maximum sensitivity of a NW- 
FET for dielectric thickness measurements or for proximity 
sensing. 


nesses much smaller than the distance nanowire-gate be¬ 
cause of the dielectric screening. We further discuss this 
property for sensing purposes. Thus, we analyze the ca¬ 
pacitance change with respect to dielectric thickness as 
an indicator of liquid height in microfluidic systems. The 
capacitance curves have an S-like shape with three re¬ 
gions. The first region is almost linear and might be 
optimum for sensing if the ratio between gate-nanowire 
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distance and the nanowire radius is about 10. The second 
region which is of the size of two diameters around the 
center of the nanowire has the largest slope and can be 
used for sensors of proximity. The third region is rather 
unimportant since the capacitance in the region saturates 
slowly to its asymptotic value given by the embedding of 
the nanowire in a dielectric of infinite thickness. 
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